Canted phase in double quantum dots 
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We perform a Hartree-Fock calculation in order to describe the ground state of a vertical double 
quantum dot in the absence of magnetic fields parallel to the growth direction. Intra- and interdot 
exchange interactions determine the singlet or triplet character of the system as the tunneling is 
tuned. At finite Zeeman splittings due to in-plane magnetic fields, we observe the continuous quan- 
tum phase transition from ferromagnetic to symmetric phase through a canted antiferromagnetic 
state. The latter is obtained even at zero Zeeman energy for an odd electron number. 
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Interaction in two-dimensional (2D) electron gases 
leads to new quantum phases when more degrees of free- 
dom (external fields, spin and layer indices) are supplied 
to the system. In bilayer v = 2 quantum Hall (QH) 
structures theoretical calculations have predicted 0] and 
experimental evidence has confirmed ||^ the existence of 
a particularly exotic canted antiferromagnetic (C) phase 
which continuously connect the naively-expected ferro- 
magnetic (F) and paramagnetic (P) ground states (GSs) 
as the layer separation is tuned. Interesting predictions 
regarding C states in few-electron double quantum dots 
(DQDs) in the QH regime have also been reported 
Thus high magnetic fields seem an unavoidable condi- 
tion to observe this quantum transition given that for 
vanishingly small magnetic fields F and C states in bi- 
layer systems would cost a good deal of kinetic energy. 
Our goal is to challenge this idea by lowering the sys- 
tem dimensionality and benefitting from the atomic-like 
spectrum of a semiconductor quantum dot In DQDs 
(termed artificial molecules as well) Coulomb blockade 
effects §, magnetization |^ , and the formation of a delo- 
calized molecular GS are some of the exciting features 
observed. From the theoretical viewpoint exact diagonal- 
ization methods ^ , Hubbard-based models (|] , and spin 
density functional theories have been developed to 
show the presence of magic-number, molecule-type, and 
Hund's-rule-violating states in vertically coupled dots. 

Here we present a Hartree-Fock (HF) theory for ad- 
dressing many-body effects in two vertically-coupled 
parabolic quantum dots separated by a distance a with a 
total electron number N . We study this system in the ab- 
sence of magnetic fields perpendicular to the dots. Still, 
in order to add spin symmetry breakings we allow for a 
parallel magnetic field whose coupling to the electronic 
orbital motion is neglected (a is assumed to be much 
smaller than the corresponding magnetic length). We 
are interested in quantum dots whose atomic-like charac- 
ter results in half-filled shells formed by quasi-degenerate 
eigenstates, thus having large spin expected values act- 
ing as effective magnets. Our main findings are: (i) 
the existence of a robust C phase (envisaged as tilted 



spin vectors) at finite Zeeman energies for even values 
of N linking the F (fully spin-polarized or, equivalently, 
triplet) and P (fully isospin-polarized or singlet) GSs via 
a second-order phase transition; (ii) the persistence of C 
states for N odd even in the absence of Zeeman gaps; and 
(iii) the overture of a simple model which qualitatively 
explains our results, allowing to deal with more complex 
quantum dot systems. 

The electron spatial coordinates are denoted by f = 
(x, p) where p = {y, z) and x is the growth direction. 
The wave function of the i-th electron may be expanded 
in terms of 2D harmonic oscillator eigenstates, (j)niip), 
n — 0,1,2, . . . being the radial quantum number and I the 
angular momentum obeying / = — n, — n + 2, . . . , n — 2, n. 
Accordingly, the Hamiltonian reads: 
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where the sums are extended to all indices, a and a are 
the spin and layer indices, and e„ — huj{n + 1) (huj is 
the confinement strength). The second and third terms 
describe the coupling of the isospin and the spin of the 
system with external perturbations, namely the tunnel- 
ing t and the Zeeman splitting Az- The isospin points 
along +{—)z when the electron is at the top (bottom) 
layer. Notice that the tunneling term only switches layer 
indices, thus conserving the rest of quantum numbers. 
Az = gpsB, where g is the Lande factor, fiB the Bohr 
magneton, and B the applied magnetic field in the z di- 
rection. V is the matrix element of the Coulomb potential 
V{\r-r'\). 

Because we seek to identify spin and particle-number 
broken symmetries which are refiected in the interdot 
coherence, we restrict our study to radial-symmetry- 
conserving solutions within an isolated dot despite 
spontaneously radial symmetry breaking states might 
take place The resulting states ip can be la- 

beled with individual angular momenta: ipi- (r) = 
J2n a a ^maia, (f>nih{p} fa,{x) Xa,- d are the coefl[icients 
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of the expansion to be self-consistently calculated, /„ (x) 
is the vertical component of the wave function, and Xo- 
is the spin part. Hence the Hamiltonian is numerically 
diagonalized in separate /-subspaces (though the matrix 
elements do depend on the total system configuration). 

We consider that only the lowest-energy states in the 
vertical structure are populated, and approximate the 
form of f{x) as follows: fa=T{x) = \/S{x) {fa=B{x) — 
\/S{x — a)) for the top (bottom) layer. More precise ex- 
pressions for /(x) would involve accurate form factors 
entering the final results but qualitatively all would 
remain the same. Then taking a close look at the x part 
of V, 

V<x fdx f dx'f^^{x)U^ix')Vi\r-r'\)UJx')U,Jx), 



we observe that the only terms different from zero are 
either ai = aj = = ctm {intradot interaction) or 
^ aj = ak {inter dot interaction) since crossed 
terms (e.g., ai — ak ^ aj — am) would be null because 
top and bottom wave functions have zero overlap. 

In the HF approach the electron-electron interaction 
part of the Hamiltonian of Eq. (|l|) can be arranged in 
two parts: the Hartree operator , 



nirimlilmf^iOi 



(2) 



where 



UjUkljlkCjOl' 



(the indices of V are omitted for the sake of simplifying 
notation) and the exchange part i?^: 

= ^ ■^ninfciiii.(Ti(TjaQ'cJj .;.^.q,C„j.;j.(j^q' (4) 



where 



/ J \ rijljCTja' 



(5) 



Throughout our calculation we make ^c[. q^^ oc 
according to the aforementioned radial symmetry ap- 
proximation. In Eq. (|^) a' = a gives rise to intradot 
Hartree interaction. The effect of this term is to make 
electrons repel from each other within the parabolic well. 
The total energy of the dot is thus augmented. Inter- 
dot Hartree interaction is naturally included for a' ^ a, 
notwithstanding it docs not have a strong influence in 
the final magnetic configurations, for it merely involves 
a rigid shift of the energy levels. However, we keep it 
for having the same number of electrons within each dot 
when N is even and for obtaining a more reliable value of 



the total energy of the system. Intradot exchange inter- 
action favors spin alignment within each dot as expected. 
For a' ^ a we are left with interdot exchange interaction. 
We stress that although the barrier separating the double 
well in the x direction is wide enough and consequently 
the vertical parts of -0 have negligible overlap (in our 
case the overlap is strictly zero owing to the Dirac delta 
functions) the interdot exchange interaction cannot be 
disregarded as it plays a crucial role in the final DQD 
magnetic order. In fact, it is the competition between 
the intradot exchange part and the interdot exchange 
interaction plus the tunneling term and the Zeeman en- 
ergy what drives the system from a GS dominated by in- 
tradot contributions (large values of the interdot distance 
or small tunneling parameter) to a state in which interdot 
effects prevail (small values of a or large t). In between 
non-trivial quantum phases can occur. Our allowance 
of significant non-zero order parameters, (^c|,j,C|b^ 7^ 0, 

^c|rpC|T^ 7^ 0, etc., eventually leads to spontaneous spin 
symmetry breakings, spin rotations, canted phases and 
the fact that the particle number at each dot is not a 
good quantum number. 

At this point a small digression about the trust- 
worthiness of the HF model is needed. Previous 
works have studied in detail |Q the differences be- 
tween (un)restricted HF theories and exact diagonaliza- 
tion methods in quantum dots. It seems clear that a 
large number of electrons would result in a negligible 
amount of quantum fluctuations which are unlikely to de- 
stroy mean-field-based predictions. Furthermore, we can 
take advantage of large values of {Sz) (highest half-filled 
shells) to ease the appeareance of C phases. Incidentally 
the existence of lower closed shells (hereafter designated 
as core) is a crucial difference between QH systems and 
DQDs in the absence of magnetic fields. In the former 
case, only the lowest Landau level is occupied and the 
kinetic energy plays a minor role. In the latter, the dot 
fills its levels following an aufbau rule, thereby closing 
shells as N is increased. In principle we could fail to take 
into account the core states and utilize them simply to 
assess the highest-lying eigenstates before coupling the 
dots (then reducing the complexity of the calculation). 
However, this would imply spin rotations applied solely 
to the last half-filled shell, giving rise to a cost of energy 
which we are unable to evaluate a priori (as it is a self- 
consistent quantity) . The total number of electrons must 
then enter the numerical procedure. 

Now, large values of N tend to contract the (renormal- 
ized) energy level interspacing in order to build a semi- 
classical radial density. This may involve the drop of 
valence electrons below the core levels and a subsequent 
reduction of {Sz)- A more favourable situation can be ac- 
complished in part by enhancing the confinement. Hence 
one should reach a compromise between these competing 
factors. 
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FIG. 1. (a) Total energy of a DQD with N = 32, 
hu} = 30 meV, a = Zo/2, and Az = 0.017?J.tJ. A dashed 
(dot-dashed) line shows the behaviour for the F (P) state 
provided spontaneous spin symmetry breaking had not been 
taken place, (b) Expectation values of the total spin for each 
dot. z (x) components are drawn in full (dashed) lines. 

The expansion of ■0 is enlarged enough, in such a way 
that the highest n state contributes less than 0.01% to a 
typical density. We present data for iV = 32 (though sim- 
ilar results are found for N — 18), setting hu! = 30 meV 
and a = lo/2, Iq = ^/fij {ra*Lij) being the harmonic os- 
cillator typical length (to* is the GaAs effective mass). 
As a result, the single-particle value of the Sz projection 
onto the layer a is (Sza) — 2. 

Fig. |l](a) depicts the total energy of the system, E — 
(7i), as a function of the tunneling parameter. At low 
t the GS is ferromagnetic, the intradot interaction being 
more important than the interdot one plus the tunnel- 
ing term. Because this is a fully spin-polarized state, 
its energy does not depend on t and remains constant 
until the system undergoes a continuous quantum phase 
transition to the C phase. In this case the system low- 
ers its energy by increasing the tunneling contribution. 
This favors the formation of singlets as well as the in- 
crease of interdot coherence (see below). To make the 
loss of Zeeman energy as small as possible the spin con- 
figuration is then canted. By further enhancement of t 
the P phase is achieved and a linear dependence of E on 
the tunneling is obtained. The C phase is thus a linear 
combination of the wave functions associated with F and 
P states. In Fig. |l|(b) we have plotted (S'za) for both 
dots. Its maximum is reached when the DQD is F, then 
it is depressed when entering the C phase. The conver- 
sion to the singlet state ((^z) = 0) is continuous and 
{Sxt) = — (Sxb)- This is the key feature of the appear- 
ance of a C phase- total spin components in the plane 



perpendicular to the field which are antiferromagnetically 
correlated. These conclusions are reinforced by examin- 
ing the interdot coherence of the top (bottom) quantum 



dot: A„ 



4| that all the Ao-o-'a components are zero in the 
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triplet phase (see Fig. |^(a)). As t increases, the system 
acquires interdot coherence until it is completely coherent 



in the symmetric state (for which A 



and 



A|fQ = ^iia 7^ 0). As C phases involve a spin symme- 
try breaking, — —^iia- Moreover, since (Ix) yields 
half the difference between the number of electrons in 
symmetric states and those in antisymmetric states, (Ix) 
is zero in the F case and reaches its maximal value in the 
P case. As Fig. ||(b) shows, the C phase develops inter- 
mediate quantities. The z isospin component would be 
different from zero provided there is more charge in one 
of the wells, e.g., by applying an electric bias; but this is 
not the present case. 
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FIG. 2. (a) Up-up and down-down spin interdot coher- 
ence (full lines) for the top layer. Up-down (dashed line) 
and down-up (dot-dashed line) spin interdot coherence is also 
shown. They are different from zero only in the C phase, (b) 
X component of the total isospin. 

In Fig. ^(a) we draw the entire phase diagram which 
characterizes the distinct GSs as a function of the Zee- 
man energy and the tunneling. For large Az and small t 
the DQD is in the spin-polarized phase. In the opposite 
limit, the singlet state energy is lower. The C phase lies 
between them. In the case of Az = we obtain a purely 
antiferromagnetic or Neel (N) GS with the spins pointing 
across from each other. 

A more striking feature is observed when a hole is in- 
troduced into the system. For N odd the highest-lying 
shells are not closed and the remaining hole is shared 
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by the two dots. From Fig. ||(b) we see that the region 
covered by the C phase is reduced at large t because the 
system takes advantage more easily of the possibility of 
tunneling by forming singlets. Likewise the F state is 
more energetically favoured at low transmissions. There 
is a range in the tunneling parameter at Az = where 
the lack of charge spontaneously induces ferromagnctism. 
But now the C state is extended even for Az = since 
a N phase cannot exist for an odd number of electrons. 

A simple model may be aimed to shed light on this 
phenomena. When the highest shell is occupied with an 
even number of electrons, it is reasonable that a Heisen- 
berg term accounts for the antiferromagnetic phase. In 
addition, the total energy must include a contribution 
stemming from the Zeeman energy which favours a par- 
allel spin alignment. Close to the F phase we propose the 
following energy functional: 



Et,Azi0) = -Jit)ST ■ Sb - 9^iBB ■ (St + Si 
= - J{t)S'^ cos e - 2AzS cos (61/2) . 



(6) 



Here we consider the total spins as classical vectorial en- 



^7 



= S, 



titles centered at each dot and assume 

being the angle spanned by both vectors. J(< 0) is a 
parameter fitted from the dependence of the total energy 
on t at Az = (one curve analogous to Fig. |l|(a)). As 
a result, J{t) is roughly quadratic for small t. Then the 
critical line marking the transition to the canted phase 
from the F state is achieved by setting E"{9 = 0) = 
which yields: Az{t) = —AJ(t)S/2. This gives the piece 
of parabola showed in Fig. |^(a). The curve crosses the 
abscissa axis at t = proving that no F state can ex- 
ist with Az = at finite t and the DQD takes on a 
N phase. When an electron is removed from the DQD 
the remaining hole prefers to keep its spin parallel to the 
rest while hopping from dot to dot. Therefore the sys- 
tem is sensitive to the particular spin orientation of the 
hole and Eq. (^ must include a term accounting for this 
fact: E E — t cos (0/2) (this is as if we had done an 
unitary rotation and kept only the diagonal terms) re- 

suhing in Az{t) = ^^^'^f^^' ■ Unhke the TV = 32 case, 
for A'^ = 31 = is a minimum of E at small values of 
t (in the interval of physical interest, i.e., [0, tt]) and the 
system remains spin-polarized (see Fig. ^(b)). Larger t 
means that — Q corresponds to a relative maximum of 
E{6) and one minimum at ^ shows up, fulfilling that 
the C phase is now the lowest energy state. Despite the 
simplicity of the model, the curves agree remarkably well 
with the self-consistent numerical solutions. 

In summary, our analysis of the GS of a vertical DQD 
based upon a mean-field framework predicts the existence 
of a canted phase for intermediate tunneling and not too 
high Zeeman energies. For a sufficiently high even elec- 
tron number (for which quantum correlation effects are 
not expected to qualitatively alter the conclusions) the C 



phase continuously connects F and P states as the tun- 
neling parameter is varied. When a hole is created in 
half-filled shells the kinetic energy of the remaining elec- 
tron promotes the F phase at small t and the C phase 
arises even at zero (arbitrarily small) Zeeman splitting. 
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FIG. 3. (a) Phase diagram for 32 electrons. Full lines cor- 
respond to the numerical calculation. The dashed line is ob- 
tained from a simple model. The infinite slope of the C-P 
line seems to be correct in the thermodynamic limit (b) 
Same as (a) for A''=31. 
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